Emission-Line Galaxies from the Hubble Space Telescope Probing 
Evolution and Reionization Spectroscopically (PEARS) Grism 
Survey. II: The Complete Sample 



Nor Pirzkal 1 , Barry Rothberg 1 ' 2 ' 3 
Chun Ly 1,4 , Sangeeta Malhotra 5 , James E. Rhoads 5 , Norman A. Grogin 1 , Tomas Dahlen 1 , 
Gerhardt R. Meurer 6 , Jeremy R. Walsh 7 , Nimish P. Hathi 8 , Seth H. Cohen 5 , Andrea 
Bellini 1 , Benne W. Holwerda 9 , Amber N. Straughn 10 , Matthew Mechtley 5 

I 

^ ABSTRACT 

00 

(N 

i— i We present a full analysis of the Probing Evolution And Reionization Spectro- 

scopically (PEARS) slitess grism spectroscopic data obtained with the Advanced 

^ Camera for Surveys on HST. PEARS covers fields within both the Great Obser- 

vatories Origins Deep Survey (GOODS) North and South fields, making it ideal 
as a random survey of galaxies, as well as the availability of a wide variety of an- 
cillary observations to support the spectroscopic results. Using the PEARS data 
we are able to identify star forming galaxies within the redshift volume < z 
< 1.5. Star forming regions in the PEARS survey are pinpointed independently 
of the host galaxy. This method allows us to detect the presence of multiple 
emission line regions (ELRs) within a single galaxy. 1162 Ha, [OIII] and/or [Oil] 
emission lines have been identified in the PEARS sample of ^ 906 galaxies down 

m 

in 

od \ 

Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD, 21210, US A 

2 George Mason University, Department of Physics & Astronomy, MS 3F3, 4400 University Drive, Fairfax, 
VA 22030, USA 



3 Leibniz-Institut fiir Astrophysik Potsdam (AIP), An der Sternwarte 16, 14482, Potsdam, Germany 



Giacconi Fellow 



5 School of Earth And Space Exploration, Arizona State University, Tempe, AZ, 85287-1404, USA 

international Centre for Radio Astronomy Research, The University of Western Australia, 35 StirlingHigh- 
way, Crawley, WA 6009, Australia 

7 European Southern Observatory, Karl-Schwarzschild-Strasse 2, D-85748 Garching, Germany 

8 Carnegie Observatories, 813 Santa Barbara Street, Pasadena, CA 91101, USA 

9 ESA Fellow, ESA-ESTEC, Keplerlaan 1, 2200 AG, Noordwijk, the Netherlands 

10 Astrophysics Science Division, Goddard Space Flight Center, Code 665, Greenbelt, MD 20771, USA 



- 2 - 



to a limiting flux of ^ 10 -18 erg/s/cm 2 . The ELRs have also been compared to 
the properties of the host galaxy, including morphology, luminosity, and mass. 
From this analysis we find three key results: 1) The computed line luminosities 
show evidence of a flattening in the luminosity function with increasing redshift; 
2) The star forming systems show evidence of disturbed morphologies, with star 
formation occurring predominantly within one effective (half-light) radius. How- 
ever, the morphologies show no correlation with host stellar mass; and 3) The 
number density of star forming galaxies with M* > 10 9 M e decreases by an order 
of magnitude at z < 0.5 relative to the number at0.5<z<0.9in support of 
the argument for galaxy downsizing. 

Subject headings: line: identification, catalogs, galaxies: evolution, galaxies: lumi- 
nosity function, mass function 



1. Introduction 



Emission line galaxies (ELGs) are systems selected by the presence of strong line emis- 
sions (e.g. Ly-a [Oil], [OIII], H/3, and Ha), usually detected using narrow band or grism 
surveys. The strong emission lines in these galaxies trace recent star formation activity, in 
contrast to the star formation history and properties of the global stellar populations that 
can be discerned using broad band observations. A significant amount of the light originat- 
ing from ELGs is contained in regions producing strong emission lines which in turns makes 
these objects easily identifiable. The emission lines of ELGs also provide a convenient way 
to determine the redshift of these objects. Since ELGs are selected on the basis of strong 
emission lines rather than continuum emission, selecting ELGs allows one to probe down to 
lower luminosity (and thus lower mass) galaxies compared to broad band surveys which tend 
to be incomplete at or above rrf galaxies at more distant epochs. Assuming that ELGs are 
spatially distributed in a fashion similar to other galaxies, they provide a powerful tool for 
tracing the star formation history of the Universe. 



The epoch < z < 1.5 discussed in this paper is important because star formation ac- 



tivity in galaxies has been observed to increase significantly as redshift increases (e.g. Madau 



et al. 


1998; 


Hopkins 


2004) 



as to whether the star formation density relation flattens or decreases, the initial increase in 
star formation implies that, at low z, some mechanism (s) must have occurred which quickly 
quenched star formation. If this was not the case, massive ellipticals today would still be 
observed to be forming many stars. There is also evidence that suggest that the inter stellar 
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medium, star formation rates and gas fractions differ between local and distant galaxies. 
Studying galaxy evolution at these redshifts therefore does not only require the ability to 
measure the star formation history of these objects, but also depends on our ability to prop- 
erly sample galaxies over a wide range of masses to alleviate as many biases as possible. 
ELGs are ideal for such work. As noted above, these objects are easily detected in surveys 
and they are efficient for probing to lower stellar masses as a function of telescope time 
required. The wavelength range of the ACS grism used for PEARS makes it possible to 
identify the strong rest-frame emission lines that are well known to be a sign of vigorous star 
formation (e.g. Ha, [OIII] and [Oil]) out to z ^ 1.5. In fact, in this paper, examining Ha, 
[OIII] and [Oil] emitters separately allows us to look at properties of star forming galaxies 
in increasing redshift ranges. When plotted separately, these three emission lines represent 
proxies for the redshift bins of < z < 0.5, 0.1 < z < 0.9, and 0.5 < z < 1.5, respectively. 



Identifying ELGs has traditionally been done using narrow band photometric filters. 
This technique has even been successfully applied to very high redshifts to detect Ly-a 
emitters ( Rhoads et al.||2001 ). However, while narrow band surveys can efficiently cover large 
fields of view down to relatively faint magnitudes, they are typically limited to very small and 
discrete redshift ranges. This can be partially alleviated using multiple narrow band filters 
(e.g. Subaru Deep Field, Ly et al.|2007 ) but continuous redshift coverage remains intrinsically 
limited in narrow band surveys. The Probing Evolution And Reionization Spectroscopically 
(PEARS) slitless grism spectroscopic survey provides an unprecedented opportunity to study 
ELGs in a way that cannot be achieved from any ground based observations. PEARS allows 
us to bypass the difficulties inherent in narrow band filter surveys (as noted above) and 
the limitations imposed by varying sky brightness and atmospheric emission lines which can 
limit ground-based grism surveys, and identify ELGs based solely on the direct detection 
of emission lines in dispersed slitless spectra. As previously shown (Pirzkal et al. 2006 



Straughn et al. 2008, 2009), this approach allows us to detect emission lines in very faint 



host galaxies, particularly sub-m* galaxies, over a very large and continuous redshift range 
and is particularly sensitive to high EW emission lines. Particular to the PEARS survey, the 
use multiple position angles on the sky allows us to identify emission lines using independent 
observations as well as to pinpoint the exact physical location of the ELRs within each ELG. 
Crucially, and in addition to this, the PEARS survey was designed to overlap with both the 
GOODS-N and GOODS-S fields, and there is therefore a substantial amount of ancillary 
data available, including very deep, high resolution broad band imaging ranging from the 
UV to the infrared bands. 



As we noted above, the redshift range (0 < z < 1.5) probed by PEARS is a critical tran- 
sition epoch, both in terms of star formation histories and morphological evolution. On one 
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hand, the PEARS grism slitless observations make it possible to efficiently identify emission 
lines, identify the corresponding ELRs and host ELGs. On the other hand, the GOODS 
ancillary data allow us to examine the morphology and physical characteristics of the ELGs. 
This powerful combination of data gives us an opportunity to examine the evolution of ELGs 
over a long period of time and over a much wider mass range than has been previously probed. 



This paper is organized as follow: Section [2] briefly summarizes the PEARS observations 
(HST Proposal 10530, P.I. Malhotra). Section [3] describes the data reduction and analysis 
of the sample, including detection, extraction and identification of emission lines, as well as 
completeness tests. Section 4.1 presents the PEARS [Oil], [OIII] and Ha line luminosity 



functions and their redshift evolution. Finally, Section |4.2| compares the properties of the 
PEARS host galaxies, such as morphology and luminosity, with the star formation properties 
discerned from the PEARS emission lines. All calculations in this paper assume H = 
70 km s _1 Mpc -1 and Qm = 0.3, Q\ = 0.7. All magnitudes are given in the AB system. 



2. Observations 

The PEARS observations were obtained as part of a large Hubble Space Telescope (HST) 
proposal (200 orbits, Proposal 10530; P.I.: Malhotra). The program used the Advanced 
Camera for Surveys (ACS) Wide Field Camera (WFC) in conjunction with the G800L grism 
filter. The G800L has a resolution of R ~ 69-131 and provides wavelength coverage of 0.55- 
1.05 /im across the entire ACS/ WFC field of view. A total of nine fields 11.65 arcmin 2 
for each field) were observed for ^40000 s (20 orbits) each, split evenly between observations 
taken at different position angles (P. A.) on the sky (typically 3 per pointing). Multiple 
P.A.s are important for identifying and masking contamination from other sources in the 
field and removing spurious pixels (e.g. cosmic rays, bad pixels, etc). Four PEARS fields 
are within the GOODS-N field ( |Giavalisco et al.||2004| ). Five PEARS fields are within the 



GOODS-S field, with one PEARS field re-observing the GRAPES/HUDF field flPirzkal et 
aLl |2004[ ). The combined areas of the PEARS-N and PEARS-S are 50.24 and 68.84 arcmin 2 , 
respectively. The fields and their location within the GOODS fields are shown in Figure [TJ 
Table [T] lists the PEARS fields positions and total exposure times. 
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3. Data Reduction and Analysis 
3.1. Detection of Emission Lines 



PEARS emission lines were initially detected directly from combined high signal to noise 
ACS grism slitless spectroscopic images. The method does not rely on and is independent 
of any imaging data or object catalog. The basic method used to identify emission lines 
in the PEARS data was described in Straughn et al. (2008, 2009). However, at that time, 
only the PEARS-S data were used. We now present the full PEARS dataset (comprising 
of PEARS-N and PEARS-S) which covers an area on the sky that is approximately twice 
as large. We have also employed a new, slightly refined version of our line identification 
pipeline and are able to reach down to a slightly lower flux level than before. The detection 
algorithm, dubbed "PEARS-2D," is very briefly summarized below and was applied to each 
of the individual 9 PEARS fields listed in Table [U 

1) All grism exposures obtained at the same P.A. on the sky are combined using the PYRAF 
task MULTIDRIZZLE (Koekemoer 2002). This produces a high signal to noise image that is 
free of cosmic-ray and detector artifacts. This image is then smoothed using a 13 x 3 median 
to produce a smooth high signal-to-noise grism image of the field. 

2) The smoothed image is subtracted from the combined Multidrizzle image of the field. 
This step essentially removes the underlying continuum emission from the dispersed spectra. 
The resulting continuum subtracted image is what we use to detect emission lines. 

3) Using SExtractor (Bertin & Arnouts 1996) and the continuum subtracted image, we 
identify point sources, e.g. emission lines, noise spikes, and detector artifacts. This step 
generates a list of emission line candidates for each P.A. 

4) Steps 1 though 3 are repeated separately for each available P.A. resulting in 3 or 4 (i.e. 
the number of available P.A.s) emission line candidates lists. 

5) Using a detailed knowledge of the instrument distortions and of the dispersion relation 
of the G800L ACS/WFC grism, and combining the data obtained using multiple P.A.s, we 
determine the location on the sky where each individual emission line originates (essentially 
where the different traces cross when projected onto the sky) as well as determine the ob- 
served wavelength of each emission line. We derive a unique wavelength for each line in each 
of the available P.A. (See Figure 2 of |Straughn et al.||2008|). 



The PEARS-2D method was applied separately to each PEARS field listed in Table [T] 
and a separate list of ELR candidates was therefore generated for each of the nine PEARS 
fields, ignoring for now the fact that some of these fields overlap slightly. The number of 
ELR candidates that was generated was controlled mainly by the detection threshold used to 
detect emission lines with SExtractor and by how much tolerance we allowed in the inferred 
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ELR positions and observed wavelengths of the emission lines. We adopted a detection of 
l.la per P. A., which for a a line detected independently using 3 different P.A.s corresponds 
approximately to a 2a detection. We also allowed for a 3 pixels tolerance in the physical 
location of an ELR (accounting for an imperfect knowledge of the trace of the ACS grism) 
and a 100A tolerance in the wavelength we derive for an emission line (again accounting 
for an imperfect knowledge of the wavelength solution of the ACS grism). While the use of 
multiple P.A.s was quite effective at filtering spurious detections, our tolerance limits and 
aggressive detection threshold resulted in the detection of fainter emission lines, at the cost 
of some false positive detections. 

Using the PEARS-2D method allowed us to generate a list of ELR candidates for each 
PEARS field that did not rely at all on any pre-generated object catalogs or pre-selection of 
target galaxies. We stress that a candidate ELR did not require the detection of a host galaxy 
in the field. PEARS-2D, with its multiple P.A.s strategy, has three immediate advantages 
over other methods that rely on observations taken at a single P. A.: First, we can expect to 
detect extremely large EWs that would not be identified through more traditional techniques. 
Second, we derive accurate locations of ELRs without assuming that the source is at the 
center of the host galaxy and we can therefore identify multiple ELRs within a single galaxy. 
Finally, the PEARS-2D line wavelength calibration is significantly more accurate. Normally, 
the wavelength reference point is tied to the location of the host galaxy (determined using 
a direct image taken in conjunction with the grism observations). However, since ELRs can 
be several half-light radii away (amounting to a non-trivial number of pixels) the wavelength 
solution of the ELR is affected by this distance from the center of the host galaxy. Every error 
of one pixel in the assumed position of the emission line feature results in a 40A systematic 
error in wavelength calibration. For large galaxies with multiple ELRs this can lead to errors 
on the order of several hundred A. Using the PEARS-2D, this error is avoided because the 
wavelength solution of the ELRs is determined independent of any information about the 
host galaxy. 



3.2. Extraction and Verification of Spectra 



We extracted spectra of the ELRs candidates identified in Section pO using the PEARS 



grism data reduction pipeline (e.g. Pirzkal et al. 2009). However, our pipeline was modified 
so that we used the positions of the ELRs candidates (i.e. the source of the emission features 
detected in our grism exposures) rather than a catalog of objects derived from the GOODS 
broad-band images. The extraction and calibration of the spectra were performed with 
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aXe (Pirzkal et al. 2001; Kiimmel et al. 2009) using optimally weighted extractions and 
an extraction width of 3x the measured emission regions sizes. Individual spectra were 
extracted from single ACS grism exposures and spectra taken at the same P.A. were then 
combined using aXeDrizzle. We thus obtained between 3 to 4 independent spectra for each of 
our ELR candidate (one per available P.A.). Wavelength and flux calibration were performed 
by aXe using the STScI provided calibration files for the ACS G800L grism. Emission lines 
fluxes were measured using gaussian fitting in each extracted spectra, resulting in 3 to 4 
measurements. We then computed the average line flux and associated standard deviation 
value for each of our emission line candidates. These line fluxes were corrected for Galactic 
Reddening assuming values from Cardelli et al. ( |1989 ). The corrections were negligible with 
(E(B-V)=0.012 for PEARS-N fields and E(B-V)=0.0078 for PEARS-S fields. 

Our list of ELR candidates reached down to very low flux levels and contained many 
false positives, we therefore had to vet each ELR and its associated emission line candidates. 
While we attempted to apply automatic techniques to accept or reject ELR and emission 
line candidates, we found it useful to manually vet spectra by eye. Authors of this paper 
manually examined and graded emission lines on a scale ranging from (very poor) to 5 
(very high). Each emission lines were graded a minimum of 3 times and by at least two 
different persons and the average of this grade was adopted. A final grade of 2 was found to 
correspond to a marginal w 2 sigma detection of spectra obtained in at least two separate 
P.A.s. 

We originally identified a total of 3705 emission line candidates using the PEARS-2D method 
described in Section [3Tj The visual quality assessment yielded a sample of 1162 marginally 
detected emission lines (529 in PEARS-N and 633 in PEARS-S). As multiple emission lines 
were sometimes detected in an ELR, the final number of ELRs was 985 (451 in PEARS-N 
and 535 in PEARS-S). 



3.3. Emission Line & Host Galaxy Identification 

Determining the nature of emission lines in the extracted spectra of a given ELR fell 
into two distinct categories: either multiple emission lines were detected or only a single 
emission line was detected. 

The resolution of the G800L grism is 40 A pixel -1 and is too low to resolve and identify 
close emission line pairs (e.g. [OIII] and H/3). However, in the redshift range of the PEARS 
survey (0 < z < 1.5), there are pairs of widely separated emission lines that allow for both 
line identification (via the ratio of the observed wavelengths of the two lines) and redshifts 
to be determined. The line pairs that we considered were: [Oil] and [OIII]; [OIII] and Ha; 
or CIII] and CIV. 
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However, most of the time spectra contained only one prominent emission line and we 
had to rely on a comparisons with photometric redshifts of the host galaxy from |Dahlen et aL 



(2010, 2012[ ). While we noted above that the PEARS ELRs were selected independently, we 



subsequently matched them with a host galaxy in the GOODS fields. We used the public ACS 
GOODS 2.0 data to generate mosaics of the GOODS fields and used SExtractor to generate 
segmentation maps and object catalogs of galaxies for these fields. In the great majority 
of cases, the RA and DEC of a PEARS ELR clearly fell within the segmentation map of a 
galaxy. In such cases the corresponding galaxy was assumed to be the host of the ELR and 
we assumed the redshift of the observed emission lines to be the known photometric redshift 
for the host galaxy. Some ELRs (~ 6% of objects with a significant emission line with a 
grade of 2.5 or larger) were found to lie beyond any galaxy segmentation maps In such cases, 
we adopted the photometric redshift of the GOODS object with the closest segmentation 
map (which not necessarily the closest object in cases of large extended galaxies). When 
comparing photometric redshifts to PEARS spectroscopic redshifts obtained using one of the 
pair of emission lines listed above, we found that in 94% of the case the redshifts we derive 



are within the 95% confidence regions given for the photometric redshifts of Dahlen et al. 



(2010, 2012) 



When only a single emission line was detected in a spectrum, we assumed that it could 
be either Ly-a, CIV, CIII], Mgll, [Oil], [Nelll], [OIII], H7, or Ha, observed at the assumed 
GOODS photometric redshift value. We identified the type of the emission line and computed 
a spectroscopic redshift for that object simply by selecting the line type that produced the 
closest match to the observed wavelength. 

We acknowledge that some of the emission lines discussed in this paper are blended 
lines, but this should have little effect on our analysis as weaker lines only weakly bias the 
fluxes and redshifts we derive. This is the case for example for [OIII] which is really two 
unresolved lines at 4959A and 5007A, Ha at 6562A which is blended with weaker [Nil] at 
6583A, while by [Nelll] we actually refers to the stronger component at 3868A). 

Some for the PEARS ELGs were found to contain more than one ELRs. Since these 
ELRs and the emission lines their spectra contained were analyzed separately, this provides 
us with a way to check the consistency of the PEARS spectroscopic redshifts. We found that 
in these cases the redshifts agreed to within z = ± 0.01, demonstrating excellent consistency. 
We briefly noted above that some of the PEARS fields overlap slightly. Some ELRs were 
therefore observed and analyzed independently as part of different PEARS fields. We found 
that 19 ELRs were observed twice (15 in PEARS-N and 4 in PEARS-S. The higher number of 
duplicate observations in PEARS-N is the result of the larger amount of overlapping between 
individual PEARS-N fields, as shown in Figure [T]). Comparing the observed emission lines 
wavelength, the emission line fluxes and finally the redshifts, we find that our errors were: 
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< 8X >= 18A; <j>= 8%; and < 5z >= 0.003, again demonstrating excellent consistency. 

For the remainder of this paper, noting that emission lines with a grade of 2 constitute a 
marginal detection, we adopt a minimum emission line grade of 2.5 or greater (corresponding 
to a line flux limit of approximately 5 x 10 -18 erg/s/cm 2 ) when deriving properties of Ha, 
[OIII] and [Oil] emitters (174, 401, and 167 lines, respectively). 



3.4. Spectroscopic vs. Photometric Equivalent Widths for ELRs at Large 

Radii 

One of the advantages of PEARS-2D is the ability to detect multiple ELRs within a 
single galaxy (see Figure [2]). However, at progressively larger radii from the galaxy center 
the contribution from the underlying continuum decreases. Since spectra are extracted at 
these large radii using small extraction windows, the measured EWs are generally larger 
than what would be derived by simply comparing the measured line fluxes to the total un- 
derlying continuum of the host galaxy. The EWs derived from narrow band imaging surveys 
generally rely on the latter method. To quantify any potential differences photometric EWs 
(EWphot) were computed using the measured PEARS line flux and the measured total host 
galaxy broad band flux. On average the spectroscopic EW (EW spec ) was ~ 3.5 x (EW p hot)- 
Histograms of the EW spec s for the Ha, [OIII], and [Oil] are plotted in Figure |3j We note 
that for the purposes of this paper an emission line is reported as a positive EW. 



3.5. Blended Emission [OIII] and H/3 Lines 

The ACS G800L grism cannot separate the [OIII] doublet (4959A, 5007A) and H/3 
(4681 A). These three lines appear blended in the PEARS spectra and H/3 can appear as 
an asymmetric feature to the stronger [OIII] line. To correct the [OIII] fluxes, each of the 
lines were fit using separate components. We assumed identical full width at half maximum 
and assumed a fixed wavelength separation for all three lines. Based on this we obtained 
estimates of the [OIII] to H/3 lines ratio for the ELRs. We found « 0.23 ± 0.25, 



which is consistent with the relative fluxes expected in star forming galaxies (Juneau et al. 

20llb. 



- 10 - 



3.6. Comparison with Previous PEARS-S Catalog 



The PEARS-S data were analyzed and published by Straughn et al. (2009) and are 



included as part of our analysis of the complete PEARS survey. While our improved analysis 
reaches down to fainter observed flux levels, which directly results in a larger number of ELRs 
being detected, our emission line list contains in excess of 90% of the emission lines listed 



in |Straughn et aL (2009). However, the number of Ha, [OIII], and [Oil] emission lines 
(with a grade of 2.5 or above) are 1.2, 1.6, and 2.7 times higher than in the |Straughn et aL 
d2009| ) catalog, respectively. This is a direct result of our ability to reliably detect fainter 
emission lines in these data, which particularly affects the intrinsically fainter [OIII] and 
[Oil] emission lines. Our use of optimal weighted extraction, as opposed to narrow box 
extraction (increasing signal-to-noise at the expense of flux completeness) as was the case in 



Straughn et al. (2009), results in more reliable flux measurements of these emission lines and 



the lines fluxes we measure are on average twice of those listed in Straughn et al. (2009). 



3.7. Completeness Simulations 

Table [i] lists the median and average line fluxes for Ha, [OIII], and [Oil] with a strong 
detection (PEARS grade of at least 2.5). In Figure [i] we show the distributions of the line 
fluxes for these three lines. The histograms are plotted as a fraction of the total for each line. 
Figure [I] demonstrates that the PEARS-2D line fluxes peak at values of ~ 10 -17 erg/s/cm 2 . 
The ACS G800L grism has an approximately flat sensitivity from ~ 6000A to 9500A, but 
our ability to recover emission lines from the two dimensional dispersed images needs to be 
carefully evaluated before we can say anything about the volume density of these sources. 
We determined the PEARS-2D detection limits using extensive end-to-end monte carlo sim- 
ulations. These steps are briefly outlined here: First, we started with the real individual 
PEARS ACS/WFC G800L exposures and artificially added a random distribution of sim- 
ulated ELRs (and simulated emission lines) to these data. We allowed for a wide ranges 
of line fluxes, host galaxies, and redshifts in these simulations and also allowed for random 
placement of multiple ELRs in each galaxy. Each of the nine PEARS field was simulated 
ten times, each time adding ELRs to 100 galaxies. 

Next, the simulated data were processed and identified using exactly the same proce- 
dures used for the real observations. The simulated spectra were then extracted using aXe 



and line fluxes measured as we described in Section |3.2[ Comparing input and output line 
lists, we found that the PEARS wavelength sensitivity is very similar to the inverse sensi- 
tivity of the ACS G800L grism with sharp cutoff below 5500A and above 9500A. This set 
the redshift limits of the PEARS survey for the Ha, [OIII] and [Oil] lines to be < z < 0.5, 
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0.1 < z < 0.9, 0.5 < z < 1.5, respectively. We also found PEARS was mainly limited 
by the intrinsic emission line flux , that ELRs with emission lines containing fluxes as low 
as 10 _18 erg/s/cm 2 could be detected and that we could detect more than 50% of emission 
lines with flux greater than ^3x 10 _17 erg/s/cm 2 . We also determined that we could detect 
emission lines with EW > 50 A. This is consistent with the observations of observed line 
fluxed shown in Figure [4} The exact fractions of lines recovered as a function of observed 
line flux for the PEARS-N, PEARS-S and PEARS-S-HUDF field (which is twice as deep as 
the other PEARS-S fields) are shown in Figure [5j A emission line grade of 2.5 corresponds 
to a line flux of approximately 5 x 10 -18 erg/s/cm 2 . 



3.8. Methods for Computing Luminosity Functions and SED Fitting 

We briefly describe the methods we adopted to determine the PEARS emission line 
luminosity functions for Ha, [OIII] and [Oil]. We also outline how we derive stellar masses 
for the ELGs containing the ELRs we identified in the previous Sections. Finally, we discuss 
the different approaches we adopted to account for internal dust corrections when computing 
the PEARS line luminosity functions. 



3.8.1. The 1/V max method 



This widely used method does not assume a shape for the luminosity function ^(L). 
However, one disadvantage is that it requires the data to be binned. The number of 
bins can impact the results. In this paper, the number of bins was determined using the 



Freedman-Diaconis rule (Freedman & Diaconis 1981), whereby the bin size is selected to be 
2 IQR(x)n~3 5 where IQR is the interquartile range of the data and m is the number of data 
points in the sample. Using the 1/V max method, the luminosity function is computed using 
the following formula: 



$(logL0 



AlogL 



E 



1 



f( Z ji L j) V j 



(1) 



where: |logL — logL^| < | AlogL; AlogL is the bin width; Vj is the maximum volume 
within which object j (observed to have a line flux of lj and to be at the redshift of zj) would 



be detected in our survey; and f(zj,Lj) is the incompleteness /(/) (see 3.7) remapped into 
absolute luminosity space L given the object's redshift Zj ) and is defined as: 
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where: Q is the solid angle of our survey (sr); V c (z) is the cosmological comoving volume at 
redshift z (in Mpc 3 ); R{z) is the normalized grism response function expressed as a function 
of object redshift; Given the redshift range {^i,^h} at which a given emission line can be 
observed by the ACS grism (i.e. observed at wavelengths 6000A < A < 9000A), the minimum 
redshift, 2^ m i n is z\ while the maximum redshift Zj 5max = min(^h, ^faint)? where Zf a i nt is the 
maximum redshift at which a line with luminosity L would remain above our minimum line 
detection threshold /threshold- Hence, Zf a i nt is the redshift corresponding to the distance of 

threshold 5 

where D L (zj) is the luminosity distance of object j. 



3.8.2. STY method 

The STY method is another commonly used method for estimating the luminosity 



function. In this one, one assumes that $(L) has the form of a Schechter function (Schechter 



1976): 



$(L)dL = $,(f rexp(-A)^ (3) 
which is characterized by the three parameters a ) $*, and L*. Following from |Sandage et] 



al. (1979), the probability of observing a given object j at redshift z with a luminosity Lj, 



is then: 



The joint likelihood can then be computed for the whole group of observed lines: 



L = UjpiLj.zj) (5) 

From this, can then determine values of a and L* that maximizes this likelihood. The 
overall normalization constant cannot be determined this way, because it cancels out in 
Equation [4| In this paper, we determined the values of a and L* by maximizing Equation [5] 
using a Monte Carlo Markov Chain approach. This allowed us to determine the most likely 
values of these two parameters as well as 95% credible intervals for these parameters. 
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was computed by integrating Equation [3] and normalizing the result so that matched the 
number of detected objects. was computed for each combination of a and L* in our 
Markov Chains to produce 95% credible intervals for the parameter 



3.8.3. Host galaxy SED fitting 



Properties of the host galaxies of the PEARS ELRs were estimated by fitting model 
Spectral Energy Distributions (SEDs) to their broadband photometric colors. The pho- 
tometric values were taken from the TFIT GOODS measurements, which include 10 and 



12 photometric bands in the GOODS-N and GOODS-S, respectively (Laidler et al. 2007 



Grogin et al. 2012). The photometry spans the observed UV ([/-band) through thermal-IR 
based on VLT, HST and Spitzer observations. The actual SED fitting was done using our 



own Monte Carlo Markov Chain SED fitting code (ttMC 2 ) flPirzkal et aLl|2012[ ) to obtain 
estimates of the stellar masses, extinction, and ages of the host galaxies. 7rMC 2 is a far more 
robust method of SED fitting than the standard x 2 algorithms because it takes into account 
a proper treatment of both error propagation and computation of confidence levels. A more 



detailed explanation of MCMC can be found in (Pirzkal et al. 2012]) and references within. 
A simple stellar population model from BC03 (iBruzual and Chariot 1120031) templates 



and a Salpeter IMF were used. While the choice of IMF and input models (e.g. BC03 or 
Maraston ( |Maraston| [2b05 ) ) can affect derived stellar masses, the effects are not the same at 



all redshifts. The detailed simulations presented in |Pirzkal et al. (2012) show that for the 
redshift range of interest here, stellar mass estimates from different models are consistent 
with each other to within a factor of a few. Other parameters obtained from SED fitting (e.g. 
extinction, metallicity and ages of the stellar population) are significantly more uncertain. 
For the purposes of this paper, we are primarily concerned with stellar mass, and to some 
extent extinction. In Figure [6] we show the distribution of stellar masses and extinctions for 
the host galaxies of the PEARS emission line sample. The mean stellar masses of the host 
ELGs are Log(mass) = 8.74 ± 0.97 M . We also estimate that the continuum extinction is 
relatively low with an average value of Av = 0.87 ± 0.88 mag. 



3.8.4- Internal Dust corrections to Luminosity Functions 

Using the current PEARS data, there was no way to directly infer the amount of in- 
ternal (to the host galaxy) dust attenuation affecting the line luminosities. We tested three 
methods for approximating dust corrections and compared the dust correction luminosity 
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functions to the ones from Ly et al. (2007). The first dust correction used attenuation values 
from the individual SED fits to the host galaxies. The second method relied on apply- 
ing an average extinction value of Ah^ = 1.0 mag (corresponding to A[on] = 1.88 mag and 
A[om] = 1-36 mag), as is commonly done in such a case (Hopkins 2004; Takahashi et al. 



|2007| ). The third method relied on a dust correction based on the somewhat more sophisti- 
cated luminosity dependent dust correction of |Hopkins et al. (2001) 

While the two first approaches are straight forward, they are a rather curve attempt at 
applying dust correction. Indeed, these approaches do not allow for the extinction values 
within ELRs to be different than the host galaxy as a whole and therefore do not allow 
for the possibility that ELRs might contain different amount of dust than the host galaxy 
as a whole. The luminosity dependent approach empirically attempts to circumvent this 
limitation. In this case, the amount of dust correction is correlated the the measured line 
luminosities. This method, unlike the first two, could in principle effect the shape of the 
luminosity functions that we measure. 

We find however that the three methods of correcting for dust have only limited in- 
fluence on the luminosity functions. Either simply shifting the luminosity functions by a 
fixed amount without affecting the slope (a) at all (as is the case when using the first two 
methods), or additionally only altering the slope (a) slightly (as is the case when applying 
a luminosity dependent dust correction). As expected, the [Oil] lines were more affected by 
dust than the [OIII] and Ha lines. However, in totem, for Ha, [OIII] and [Oil] the slopes 
varied by only ~ 0.1 on average using either of the three methods discussed above, and 
we conclude that the effect of dust correction on the slope of the luminosity functions is 
therefore negligible. 



4. Results 
4.1. The Emission Line Regions 

Jf.A.1. Star Forming Galaxy Density 

Using the uncorrected PEARS lines listed in Table [4j we can compute a space density of 
star forming galaxy (SFG), as measured by the PEARS survey, and compare it to previous 
ACS grism based surveys. We estimate the star forming galaxy density at0.3<z<1.3to 
be 4.5 x 10 -3 Mpc -3 . This is in complete agreement with previous pure parallel ACS grism 



surveys such as the one described in Drozdovsky et al. (2005). 
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4-1-2. Luminosity Functions 



We computed the luminosity functions for the Ha, [OIII], and [Oil] samples using both 
the 1/V max method and the STY methods that we described in Sections 3.8.1 and 3.8.2, The 
non dust-corrected luminosity functions, computed using the 1/V max method are shown in 
Figures 0to[9| In these fi gures, the new measurements are compared to those of (Ly et al. 



2007, shown with open triangles and also uncorrected for dust). We added one additional 



constraint to the data from Ly et al.] ( 2007[ ). The comparisons were made only with objects 
from that paper with an EWs > 50 A. 

The results from PEARS agree fairly well with earlier results, although they probe lower 
line luminosities for [OIII] and H<x The 1/V max results from PEARS-N and PEARS-S are 
plotted separately in Figures [7] to [9] The differences between these two large and independent 
fields are well within the errors. Table [5] summarizes the results from fitting the luminosity 
functions to each of the emission lines in each field separately and PEARS-N and -S together 
for both the 1/V max and STY methods. Table [5] also includes the associated 95% credible 
intervals. 

When using the 1/V max method, as previously noted in Section 3.8.1, the choose of bin 



size is important. We illustrate the effect of various bin sizes by showing (using light shaded 
circles) the luminosity functions we compute while allowing the bin sizes to vary. As one can 
visually witness, the effect of bin six has an immediate effect on the values with derive at 
a given luminosity. We show the luminosity derived using the optimal bin sizes using solid 
symbols. Error bars associated with individual points were derived using a few thousands 
bootstrapping iterations. 

Here, we briefly note that the STY method produces steeper slopes than 1/V max for 
[Oil]. Yet for [OIII] and Ha the two methods are consistent with each other. The differences 
between the two are likely due to the limited number of sources over a wide redshift range. 
This underlines the difficulties in obtaining luminosity function estimates with limited num- 
bers of sources as well as the impact of different methods for deriving those estimates. We 
further note that the results shown in Figures [7] to [9] are also generally consistent with those 



from Ly et al. (2007). 



4.1.3. Do Luminosity Functions Evolve with Redshift? 

Two advantages of the PEARS-2D sample are: that it reaches to faint line fluxes (a 
few times 10 -18 erg/s/cm 2 ); and its the wide redshift range coverage. These can be used to 
investigate whether or not the slope of the luminosity functions changes with redshift. A 
change in the slope of the luminosity function as a function of redshift is a strong indication 
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of evolution of how star formation has being occurring in the past « 9 billions years. To 
examine this issue, each of the emission lines was divided by redshift, first into two bins, 
then into three bins. We two different number of bins were used to provide for a consistency 
check of our results. A luminosity function was fit to the data in each redshift bin (see Table 
5j). This allowed us to compute slopes for different redshift ranges. These slope are plotted 



against redshift in Figure 10, The "resolution" (number of bins) in Figure 10 is limited by 
the available number of sources. Figure [10] shows a significant decrease in the value of a 
as redshift increases for Ha, [OIII], and [Oil] at both "resolutions." This is consistent with 
earlier results from Ly et al.| (2007). It suggests an evolution in ELGs and their capability to 
form stars over the last several Gyr as the proportion of high to low SFR regions increases 
as redshift increases. 



4-1-4- The Spatial Distribution of ELRs 

A major difference between PEARS-2D and other ELG studies is that we are able to 
detect the presence of multiple ELRs within a single galaxy. A breakdown of the sample 
shows that 69% of the ELGs contain a single ELR; 24% contain two ELRs; 4% contain 
three ELRs; and 3% contain four or five ELRs. Comparing the location of the ELR, as 
parametrized by their distance to the center of the host ELG normalized by the half-light 
radius of the host ELG, we can compare the distribution of single emission line regions to 



the distribution of multiple emission line regions. Figure Y\_ shows there is no indication 
of any strong differences between the two samples. In both cases, the emission line regions 
appear to be predominantly located around one half-light radius from the center of the 
ELG. We compute a Kolmogorov-Smirnof test (KS) p-value of 0.49 which indicates that 
these distributions are statistically identical. 



4-1-5. Star Formation Rates of ELGs 

One of the most important question in galaxy evolution is whether the SFR changes 
over time, and if so, what shape does SFR vs redshift have. Assuming that ELGs are 
representative of star forming galaxies in general, the depth of the PEARS-2D study and 
the large redshift range allows us to tackle this important question. The SFR was calculated 



for the [Oil] and Ha emission lines using the Kennicutt (1998) relations. For [OIII] which 



are likely to actually be blended [OIII] and H/3 the relation from Equation 5 in Drozdovsky 



et al. (2005) was adopted. The results are shown, as a function of redshift in Figure 12 



Also plotted (using a black solid line) is the SFR for emission lines with an observed flux of 
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3 x 10 17 erg/s/cm 2 . This corresponds to the PEARS-2D 80% completeness level (solid black 
solid line). This illustrates our ability to detect emission line uniformely from < z < 1.5. 

Simply computing the SFR does not provide an entirely accurate assessment of star 
formation because the mass of the galaxy affects the rate at which stars form. Instead, one 
can normalize the SFR by the mass of the galaxy (in this case estimated from the stellar 



masses computed in Section 3.8.3) to derive the specific SFR (sSFR). A histogram of the 



sSFR for the PEARS-2D ELGs is shown in Figure 13 



As the sSFR value of a galaxy can be taken, as a rough approximation, as one over 
the built up time of a galaxy, assuming constant star formation rate, the PEARS ELGs 
sSFRs implies a possible stellar mass built-up time of a few billions years. Note however 
that these sSFR estimates should be considered to be lower limits since some non detected 
star formation might be present in the PEARS ELGs and the SFR values we quote are not 
corrected for extinction. 



There has been some discussion ( Guo et aL]|2011 ) as to whether the SFR in star forming 
regions of galaxies should be correlated to the location of the star forming regions within 



the galaxy. We investigate this possible relation using the PEARS ELG sample. Figure 14 



shows a plot of the estimated SFR of each ELR in the PEARS-2D sample, separating Ha, 
[OIII] and [Oil] emitting regions, as a function of radial distance of the ELR (normalized 
to the half light radius of the galaxy). As this Figure shows, we see no indication for trend 
as a function of ELR location for either of the three types of ELGs we examine. A simple 
Pearsons linear correlation test for Ha, [OIII] and [Oil] yields values of 0.04, -0.001 and 
-0.05, respectively, indicating no statistical correlation between location of ELRs and SFR 
in those ELRs. 



Finally, Figure [15] compares Log SFR against Log M* for the ELGs in the PEARS-2D 
sample (open circles in all panels). The ELGs are plotted in four redshift bins to match the 
work of |Noeske et al. (2007). In that work Noeske et al. (2007) derived a "main-sequence" 
of star forming galaxies for field galaxies in the Extended Groth Strip, complete to Log M 
~ 10.8 (Figure 1 in that paper). The red squares shown in Figure 15 are the median values 
for the galaxies in Noeske et al. (2007) along with the ±la (dotted red line). The conclusion 
was that there exists a gradual decline in SF of most galaxies since z~l. The implication is 
that the same physics that regulates SF in local disk galaxies is occurring, which could be 
either an evolution in the gas supply or changes in the SF efficiences. Noeske et al. (2007) 
suggested that the slope of the "main-sequence" is related to the gas exhaustion of galaxies 
and is related to the age of the galaxy and SF timescale, all of which are dependent on the 
galaxy mass. The PEARS-2D sample probes galaxies to much lower masses than those in 
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Noeske et al. (2007). As such, we compare our galaxies to those in in Noeske et al. (2007) in 



Figure [T5} From this we can conclude that the relation we observe between SFR and host 
galaxy stellar mass is consistent with the star formation "main sequence" of objects with 
stellar masses of ~ 1O 1O M while extending this relation by about four orders of magnitudes 
in mass. We also potentially witness a flattening of the sSFR versus mass relation for lower 
mass objects (below « 1O 8 M ). The dashed horizontal lines in Figure 15 show our sensitivity 
limits at the lower and higher ends of the redshift ranges we show and it is clear that the 
flattening of this relation is not caused by incompleteness, especially at the higher redshift s. 



4.2. The Host Galaxies 



4.2.1. Morphologies 

The PEARS-2D galaxies comprise a remarkably robust sample to test the evolution 
of ELGs and compare their morphologies with properties such as SFRs and stellar masses. 
Unlike many morphological studies, our sample was not pre-selected by redshift or luminosi- 
ties. This sample serves as a random selection of star-forming galaxies within the PEARS 
redshift volume, which is mostly unbiased by the actual morphology of these star-forming 
galaxies. In this section, we parameterize the morphologies of the host galaxies using the 
Gini Coefficient G and M 2 q parameters (Lotz 2004). The G and M 2 q parameters can be 



thought of as proxies for dumpiness and concentration coefficients and have been shown to 
be a good way to distinguish between "normal" galaxies and galaxy mergers in the local 



Universe, as demonstrated by Lotz] (2004) in the blue using Sloan (/-band, £>j, Thuan-Gunn 
g, and £>-band (all corresponding to 4300 — 4500A), as well as at ~ i?-band. Local spiral 



and elliptical galaxies follow a well defined sequence (e.g. Figure 9 in Lotz ( 2004[ )), while 
mergers have larger G and smaller M 2 o values ( Lotz||2004[ 2008, 2010). In order to compare 
the PEARS ELGs to galaxies in general, we computed the rest-frame G and M 2 o coefficients 
for both the PEARS ELGs as well as the entire GOODS catalog using the GOODS 2.0 ACS 
public data. These values were measured in all available observed wavelengths and a A ~ 
B-band (or ~ 4350 A) rest-frame value was obtained by linearly fitting these measurements. 

As a comparison, we have also included field galaxies from within the same GOODS 
fields using the GOODS 2.0 ACS public data. For these objects, photometric redshifts were 
used to derive rest frame B band values for G and M 2 o- The rest-frame values of G and 
M 2 o were computed linearly fitting the values measured in each of the available bands. The 



galaxies are all plotted in Figure [16j which is divided into three rows of three panels for 
clarity. The field galaxies from GOODS are plotted as contours, and in each panel the Ha 



[OIII] , and [Oil] ELGs are shown separately. The solid lines in Figure 16 delineates disturbed 



- 19 - 



galaxies (above the line) from "normal" galaxies (below the line) according to Lotz (2004). 
We also show the regions containing early type and late type "normal" galaxies. When 
compared to the rest of the GOODS galaxies (black contours), the PEARS ELGs clearly 
have higher G and M 2 o values and fall above the fiducial line separating quiescent galaxies 
and active galaxies (from Lotz 2004). This strongly suggests that the PEARS ELGs have 
perturbed morphologies, likely due galaxy interactions. 

Figure [16] implies that the PEARS-2D ELGs have disturbed morphologies, likely from 
some type of interaction. We also see an indication that ELGs with more than one detected 
ELR tend to have more disrupted morphologies, as shown in the three panels in the middle 
row of Figure 16, while ELGs with a single ELR tend to lay closer the the line separating 
normal galaxies and mergers (bottom three panels), as defined by Lotz (2004). However, 
there is no correlation between the Gini-M20 values and computed SFRs and stellar masses. 
To test for any correlation correlation, the Pearson Correlation coefficient (r) was used. It 
tests the degree of linear correlation between two independent data sets, r ranges in value 
from -1 to +1 (perfect negative or anti-correlation to perfect positive correlation). The most 
correlated relation we find is that of the [OIII] versus stellar mass, shown in Figure [TtJ which 
is very weak with a value of r = 0.16. All other relations show no statistical correlations. 



4-2.2. 4350A Rest-frame Luminosity of ELGs 

The underlying host galaxy luminosity may provide additional information about the 
nature of the ELGs and how they compare to other galaxies in the field. Absolute magnitudes 



at 4350A (M4350) were computed for both the ELGs and GOODS field galaxies. Figure 18 
shows a histogram distribution of M4350 for the ELGs, divided into three panels, one for 
emission line. The median M4350S as a function of the emission line are: -21.2 for [Oil]; -19.0 
for [OIII]; and -18.2 for Ha. 

One important question is whether the ELGs are representative of other galaxies within 
the same volume. As discussed in the Introduction, ELGs are incredibly useful for probing 
the evolution of the SFR not only out to more distant epochs, but down to fainter luminosities 



(and thus lower masses) than other galaxies. Figure 19 compares the luminosity function 
of the ELGs (separated by emission line which, again, we point out approximates different 
redshift bins) compared with field galaxies from GOODS, also delineated by redshift. Figure 



19 clearly shows that the maximum density of Ha, [OIII], and [Oil] occurs at much brighter 
levels than the depth of the GOODS data. While the luminosity functions of the GOODS 
field galaxies increases monotonically in a Schechter-like manner, the number of low mass line 
emitters decreases quickly. It is unlikely this trend is due to incompleteness. The PEARS-2D 
method is more sensitive to higher equivalent width lines, therefore we would expect to see 
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more emission lines with observed flux levels of 3 x 10 -17 erg/s/cm 2 when the continuum 
light from host galaxies decreases. Yet, there appears to be an absence of faint emission line 
galaxies in all three sub-samples. 

To further quantify the effect, we examined the volume densities of [OIII] emitting 
galaxies (our largest sample) at redshifts from to 0.9. This is the redshift range with 
the largest number of ELGs. We imposed a further restriction and only selected galaxies 
above the completeness limit. This group was then divided into two subsets (0 < z < 0.5 
and 0.5 < z < 0.9) to see if the trend remained. The results are plotted in Figure 20 for 



^4350 (left) and Log Stellar Mass (right). The left panel of Figure [20] confirms that there 
appears to be a relatively small number of faint galaxies with detected [OIII] emission at 
higher redshifts. Recall that these host galaxies were selected solely based on the direct and 
independent detection of [OIII] in emission and thus independently of their observed size 
and host luminosity. 

In the right panel of Figure 20 the stellar masses are compared for the two redshift 
ranges (same limits on sample selection as the left panel). The stellar mass distribution of 
galaxies with detected [OIII] emission differs significantly. At lower redshift there appears 
to be fewer massive galaxies with detected star formation. We conclude that, strongly star 
forming galaxies were on average more massive at higher redshifts. These results seem 



consistent with downsizing (Cowie et al. 1996, e.g. 



5. Conclusion 



We have presented a sample of ELGs selected independently by their emission lines 
without a priori knowledge of their host properties. The methodology used (PEARS-2D) 
is based on direct detection of emissions line from HST slitless grism spectroscopy, with the 
added bonus of being able to detect multiple ELRs within a single galaxy. This has allowed 
us to construct a sample which is effectively random, and blind to other parameters. Using 
the wealth of ancillary data, we then investigated and compared the properties of the under- 
lying hosts with the SFR histories derived from the ELRs. The key results are summarized 
below: 

1) There is evidence for evolution in the luminosity function of Ha, [OIII] and [Oil] emission 
lines. The slopes flatten as a function of redshift. 

2) The morphology of the host galaxies clearly indicates that these objects are disturbed, 
although we detect no correlation between morphology and our stellar mass estimates, or 
star formation intensity. 

3) The mass density function of [OIII] emitting galaxies at < z < 0.9 strongly decreases. 
The number density of objects with stellar masses greater that « 10 10 M undergoing strong 
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SF decreases at lower redshifts. This supports the idea of downsizing (Cowie et al. 1996 
e.g.). 

The results presented here also demonstrate the clear advantage of space-based grism spec- 
troscopy using multiple position angles. Such observations are able to probe deeper than 
similar ground-based studies. The PEARS-2D method also provides a method for detecting 
multiple ELRs and allows spatial information about SF to be derived for galaxies. Future 
work will include using the WFC3 grism mode, with observed wavelength coverage of 0.8- 
1.6/xm. This will allow us to probe to significantly higher redshifts and determine whether 
the trends reported here continue to earlier epochs. 

Acknowledgments - NP wishes to thank Kai Noeske for the stimulating discussions. This 
research made use of the OSX version of SCISOFT assembled by Dr. Nor Pirzkal and F. 
Pierfederici. 
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Table 4. Properties of significantly detected (grade>2.5) emission lines in the PEARS 

sample 
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Fig. 1 .— The location of the four PEARS-N (left) and five PEARS-S fields (right) within 
the GOODS-N and GOODS-S fields. The fields are oriented so that North is up and East 
is to the left. Each of the shown PEARS field is approximately 200" arc second wide. Note 
that the total area where PEARS fields overlap is higher in PEARS-N than in PEARS-S. 
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Fig. 2. — A sample of PEARS star forming galaxies with their identified star forming regions 
by red squares ([Oil]), green triangles ([OIII]), and blue circles (Ha). The redshift is indicated 
at the top left of each stamp image and the one arc second scale is shown at the bottom left 
of each stamp images. 
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Fig. 3. — Rest-frame spectroscopic EW in the PEARS-2D sample. The median EW spec we 
measure are 43A, 89A, and 69Afor [Oil], [OIII], and Ea emission lines, respectively. These 
lines act as proxies for the redshift ranges of < z < 0.5, 0.1 < z < 0.9, and 0.5 < z < 1.5, 
respectively. 
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Fig. 4. — Distributions of observed line fluxes in the PEARS-2D sample, uncorrected for 
completeness or dust extinction. These lines act as proxies for the redshift ranges of < z < 
0.5, 0.1 < z < 0.9, and 0.5 < z < 1.5, respectively. 



- 31 - 



i | i | i | i | i | i I 




Log(flux)(erg/s/cm 2 ) 

Fig. 5. — PEARS sensitivity to emission line flux. Based on our simulations, we estimate 
that we can reliably detect emission lines with fluxes greater than 3 x 10 -17 erg/s/cm 2 over 
the whole PEARS fields while the PEARS-S-HUDF field, which was observed twice as long 
as each of the other 8 fields, reaches line fluxes approximately 1.4 times fainter. 
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Fig. 6.— Left panel: Histogram of the PEARS emission line galaxies stellar masses 
(in Mq) as determined from SED fitting. We derive a mean stellar mass of 
Log(mass) = 8.74 ± 0.97 M . Right panel: Histogram of the PEARS emission line galaxies 
extinction (Av in mag) as determined from SED fitting. We derive a mean extinction of 
Av = 0.87 ± 0.88 mag 
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Fig. 7. — PEARS Ha luminosity function at < z < 0.5. We show the 1/V max results 
for the full PEARS, PEARS-N and PEARS-S in black circles, blue upright triangles, and 
red downward triangles, respectively. The fits to the 1/V max results are shown by the solid 
black line, blue dashed line, and red dash-dotted line, for PEARS, PEARS-N and PEARS-S. 
We see no significant differences between the PEARS-N and PEARS-S fields. We also plot 
the sample of z = 0.4] Ha emitters from Ly et al. (2007), also with no dust correction and 
excluding objects with EW < 50 A from their sample so that we can better compare our 
results and demonstrate how PEARS reaches down to fainter luminosities. 
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Fig. 8.— PEARS [OIII] luminosity function at 0.1 < z < 0.9. We show the 1/V max results 
for the full PEARS, PEARS-N and PEARS-S in black circles, blue upright triangles, and 
red downward triangles, respectively. The fits to the 1/V max results are shown by the solid 
black line, blue dashed line, and red dash-dotted line, for PEARS, PEARS-N and PEARS-S. 
We see no significant differences between the PEARS-N and PEARS-S fields. We also plot 
the sample of z — 0.6 [OIII] emitters from Ly et al. (2007), also with no dust correction 
and excluding objects with EW < 50 A from their sample so that we can better compare our 
results and demonstrate how PEARS reaches down to fainter luminosities. 
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Fig. 9. — PEARS [Oil] luminosity function at 0.9 < z < 1.5. We show the 1/V max results 
for the full PEARS, PEARS-N and PEARS-S in black circles, blue upright triangles, and 
red downward triangles, respectively. The fits to the 1/V max results are shown by the solid 
black line, blue dashed line, and red dash-dotted line, for PEARS, PEARS-N and PEARS-S. 
We detect no significant differences between the PEARS-N and PEARS-S fields. We also 
plot the sample of z — 0.9 [Oil] emitters from Ly et al. (2007), also with no dust correction 
and excluding objects with EW < 50 A from their sample so that we can better compare our 
results. 
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Fig. 10. — Slope of the PEARS [Oil], [OIII], and Ha luminosity functions as a function of 
redshift. In the left panel, each set of emission lines is split into two groups of redshift ranges 
(yielding 6 bins). This shows a clear trend in which the slope (a) decreases with increasing 
redshift. In the right panel, each set of emission lines was divided into three redshift bins 
(yielding 9 bins in total). Again, the trend remains present. We present both a two and a 
three bins solution to show that the trend we show is not strongly affected by the binning. 
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Fig. 11. — Histograms of the distance of the PEARS ELRs from the center of their host 
galaxies, measured in units of multiple Ri. We show the distribution of ELRs in ELGs 
where only one ELR was identified as well as the distribution of ELRs in ELGs where more 
than one ELR was identified. There is no evidence that these distributions are different and 
a K-S test p-value of 0.49 confirms this. 
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Fig. 12. — SFR as a function of redshift for the PEARS [Oil] [OIII] and Ha sources with a 
line grade greater than 2.5. The solid black line shows the SFR corresponding to a flux limit 
(~ 80% completeness limit) of 3 x 10 -17 erg/s/cm 2 . The dashed black line corresponds to 
an emission line flux of 1 x 10 -17 erg/s/cm 2 
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Fig. 13 — Distribution of sSFR for the PEARS host galaxies. 



- 40 - 




Fig. 14. — SFR of the PEARS emission line regions plotted as a function of their radial 
position in their host galaxy, normalized to the half light radius of the host galaxy). The 
amount of star formation appears completely uncorrelated to the location of the ELR in the 
host galaxy with Pearsons correlation coefficients of -0.001, -0.05 and 0.04 for Ha , [OIII] 
and [Oil], respectively. These three panels are proxies for the redshift bins of < z < 0.5, 
0.1 < z < 0.9, and 0.5 < z < 1.5, respectively. 
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Fig. 15. — Comparison of SFR with Stellar Mass for the ELG sample (open black circles). 
The data have been binned into 4 redshift ranges for 1:1 comparison with higher mass star 
forming galaxies from the Extended Groth Strip (Noeske et al. 2007). The dash-dotted red 
line represents the ±la of the median values. The solid black circles represents the median 
for the ELG sample. The solid black lines represents the ±la of the median values. The 
dashed horizontal line represents the 80% completion levels at the minimum and maximum 
redshifts considered in each panel. 
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Fig. 16. — Morphology of the PEARS ELGs as parametrized by the Gini and M20 coefficients 
at the rest-frame wavelength of 4350A. Top row: PEARS ELGs with [Oil], [OIII], and Ra 
ELRs (red square, green triangle and blue circle, left to right columns, respectively). Middle 
row: PEARS ELGs containing multiple ELRs. Bottom row: PEARS ELGs containing 
only one ELR. In every panel we show the rest-frame morphology of the GOODS field 
galaxies with photometric redshift estimates that are within the appropriate redshift ranges 
< z < 0.5, 0.1 < z < 0.9, and 0.5 < z < 1.5, left to right columns, respectively. The hosts 
of the emission lines we detected are nearly all above the line (shown in black) separating 
"normal" galaxies (below the line) and "merging" galaxies in the nearby Universe and is 
taken from |Lotz| ( |2004[ ) . 
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Fig. 17. — Gini coefficient values of the PEARS host galaxies versus their stellar masses, as 
estimated from SED fitting. The M20 and the Gini values are shown in the top and bottom 
row, respectively. We also show the Ha, [OIII] and [Oil] host galaxies separately in the 
left, middle and right most column respectively. There is little evidence for a strong trend 
between stellar mass and either the M20 or Gini coefficients in our PEARS emission line host 
galaxies, as indicated by Pearsons correlation coefficient values of at most w 0.16. However, 
a mild decrease in M20 and increase in the Gini coefficient as stellar mass increases can be 
seen for the [OIII] and Ha host galaxies (at redshifts of0.1<z<0.9 and < z < 0.5 
respectively) . 
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Fig. 18. — Distributions of the 4350Arest-frame absolute magnitude of the host galaxies of 
the PEARS Ha, [OIII] and [Oil] emission line (redshift ranges of < z < 0.5, 0.1 < z < 0.9, 
and 0.5 < z < 1.5, respectively). 
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Fig. 19. — Luminosity functions of the 4350Arest-frame absolute magnitude of the host 
galaxies of the PEARS Ha, [OIII] and [Oil] emission line (symbols with error bars, redshift 
ranges of < z < 0.5, 0.1 < z < 0.9, and 0.5 < z < 1.5, respectively). The luminosity 
function of GOODS field galaxies at respectively similar redshifts are shown (scaled) using 
solid lines. The density of Ha emitters (0 < z < 0.5) peaks at « —18, while the density of 
[OIII] emitters peaks at « — 19 and the density of [Oil] emitting galaxies (0.5 < z < 1.6) 
peaks at « —21. We show both the completeness corrected (filled symbols) and uncorrected 
(open symbols) density estimates. 
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Fig. 20. — Left Panel: Luminosity functions of the 4350A rest-frame absolute magnitude 
of the host GOODS galaxies where at least one [OIII] emission line region was detected 
with a corresponding SFR of 1.7 M yr _1 or greater (corresponding to our flux limit at 
the maximum observable redshift of 0.9). We plot the densities of host galaxies in the 
lower redshift range of 0.1 < z < 0.5 using green circles, and the ones within the redshift 
range of0.5<z<0.9 using red triangles. There is a strong decrease in the density of faint 
(M4350 > —19) host galaxies at higher redshifts that does not exist at lower redshifts. Right 
Panel: The corresponding host galaxy mass function for the data shown on the left panel. 
There is a strong, ten-fold decrease in the number density of galaxies with [OIII] emission 
with stellar masses greater than 10 9 M in the lower redshift bin while the density of [OIII] 
emitting galaxies remain the same for less massive galaxies. 
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